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Abstract 

We estimate an unknown qubit from the long sequence of n ran- 
dom polarization measurements of precision A. Using the standard 
Ito-stochastic equations of the aposteriori state in the continuous mea- 
surement limit we calculate the advancement of fidelity. We show 
that the standard optimum value 2/3 is achieved asymptotically for 
n ^> A 2 /96 3> 1. We append a brief derivation of novel Ito-equations 
for the estimate state. 

1 Introduction 

The standard object of quantum inference is the value a of some hermitian 
observable a of the given quantum system. The process of inference is called 
quantum measurement. One can consider the apriori quantum state p of 
the given system as an additional object of inference |]l|, M. The limitations 
as well as the optimization of state determination are in the focus of recent 
investigations ]3], |||, [| especially in the field of quantum information and 
communication ||. A completely unknown state p can not be inferred from 
a single system: the fidelity of the estimate p' will be poor. If the apriori 
state p is pure then the estimate p' must also be pure, and the simple 
bilinear expression F = trf/S'p] defines its fidelity. If we assume that the 
apriori pure p is completely random then lower and upper limits become 
analytically calculable for the average fidelity F §. For a single two-state 
system (qubit) one obtains: 

- < F < - . (1) 
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Any deliberate trial p' , when completely unrelated to p, will yield the same 
worst value 1/2. The best value can be attained in many ways. Let us, for 
instance, measure the Pauli-polarization matrix a along a single randomly 
chosen spatial direction. Let a = ±1 be the results of the projective mea- 
surement. It is then natural to identify the estimate pure state p'(cr) with 
the standard aposteriori pure state p(a) taught in textbooks: 



Trivial calculation can prove that the average fidelity over random apriori 
pure states p is 2/3. 

No quantum measurement however involved could improve on F = 2/3. 
In particular, it would make no sense to perform a second projective mea- 
surement on the given single qubit. We can, however, consider non-projective 
measurements M, |j from the beginning. A typical non-projective measure- 
ment yields less information than an ideal measurement would do. Hence it 
makes sense to combine successive non-projective measurements on a single 
system || in order to improve fidelity. In what follows, we mean non- 
projective measurements unless we say otherwise. 

The general case involving a sequence of repeated measurements is be- 
yond the capacity of analytic calculations. There is, nonetheless, an effective 
theory for long sequences. Then the measured value a, the aposteriori state 
p(c), and the state estimate p'(o~) all become time dependent and satisfy 
coupled stochastic differential equations. The 'conditional' master equation 
of the aposteriori state Q as well as its coupling to the measured value [fL0| l 
have been well-known from long ago (see also pi] ]) as the ultimate formal- 
ism of earlier continuous measurement models |13| . The equation of the 
estimate state has remained undefined and we outline its derivation in the 
Appendix. 

In Sect. 2 we discuss state estimate from a single measurement. We suc- 
ceed to express the average fidelity in terms of aposteriori states. In Sect. 3 
this result is generalized for a sequence of measurements. In Sect. 4 the con- 
ditional 'master' equation is introduced for the aposteriori state. In Sect. 5 
we calculate the progression of fidelity for long sequences of very unsharp 
measurements and we prove how fidelity will saturate to 2/3. Although 
we develop the concrete equations for two-state systems, most results can 
trivially be extended for higher dimensions N. 
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2 Fidelity from single measurements 



We approximate the exact eigenstates of a given hermitian observable a by 
approximate Gaussian projectors of precision A: 



1 



: exp 



They satisfy the completeness condition 



cr 



2A 2 



U(a)da = I 



(3) 



(4) 



and form a POVM |7|, g]. In the simplest case, the corresponding (non- 
projective) measurement of a will transform the apriori state p into the 
following aposteriori state: 



tr 



U(a)p 



(5) 



where cr is the random outcome of the measurement. It may take any real 
value with the normalized probability density 



p(a) 



tr 



U(a)p 



(6) 



The theory of (non-projective) measurements does not imply a theory 
for the estimate p'. One could mistakenly think the aposteriori state p(a) a 
reasonable estimate for the apriori state p. Unfortunately, the experimenter 
has no access to it. He/she infers the measured value a and it is, contrary 
to the projective measurement (||) , not enough to derive the aposteriori 
state. It is only sufficient to identify the approximate projector n(<r). Its 
normalized form can be a reasonable estimate: 



n(g) 

trn(a) 



(7) 



This is a mixed state. If the apriori states p are unknown pure states then 
the estimate should also be pure. To this end, the experimenter must refine 
his/her first choice (R) . The estimate will be one of the pure eigenstates of 
the mixed state estimate (^) , chosen randomly with probability equal to 
the corresponding eigenvalue. (The optimum estimate would be the most 
probable eigenstate M.) 
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In our work, we discuss pure apriori states and, accordingly, we use 
the above mentioned pure state estimates. In other words, the pure state 
estimate will be an eigenstate of II(cj), with probability proportional to the 
corresponding eigenvalue of tl(a). By construction, the average of these 
pure state estimates is identical with the mixed state estimate (|7j) . This 
has a useful consequence in fidelity calculations. The bilinearity of fidelity 
tr[p'p], valid originally between two pure states, will be preserved for the 
expected fidelity of our estimates: 

F = J tr [p'(a)p] p(a)da = E \x\pt{o)p] , (8) 

where p' is defined by (0) and E stands for stochastic expectation value. 

We benefit from the bilinearity. We are going to find a simpler expression 
for F. While we retain the notation p for the pure apriori state, we imagine 
a hypothetical apriori state fr = I/N as well, which is totally mixed. We 
apply the non-projective measurement (|5],||) to jr . This yields the simple 
relationship n(cr) = Np~ (cr)p- (<r) where p ? (cr) = iV _1 trII(<7) is the proba- 
bility distribution of the outcomes for the measurement on the hypothetical 
apriori state jr . Substituting these relationships into (Q) and inserting the 
result into @ , we obtain the following new form: 

2 



N / tr 



pV)pjj p'(a)da = NE[tT[pi r (a)p\) . (9) 

Note that the stochastic average is to be taken with the hypothetical prob- 
ability distribution p ? (cr) instead of the true p(<r). The new expression (^) 
contains the (hypothetical) aposteriori state while the old formula (|8|) con- 
tained the (true) estimate state. It pays because the aposteriori states will 
satisfy simpler equations than the estimate states, see Sect. 4 and the Ap- 
pendix. 

If we follow the example of Sect.l, we have to average the fidelity (|) 
over random pure qubit states p: 

F = i + ^E tr[pV)] 2 . (10) 

This formula of the average fidelity will be generalized for the continuous 
estimation of random pure states in Sect. 3. 

3 Fidelity from sequential measurements 

We start from the sequence 111 [a\ ),..., II ra (cr n ) of n measurements 

The measured observables &i,...,a n need not to be identical. Thus our 
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measurements may not commute. It is well-known that a sequence of mea- 
surements is formally equivalent with a single (though complicated) mea- 
surement. Applying the Eq. (H) n times repeatedly, the aposteriori state 
becomes: 

G n (a)pGl(a) 



P 



tr 



n n (cr.)p 
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The sequentially composed Kraus- operator 0, |(| reads: 

G n {a)=m\a n )...h\ l \a 1 ) , (12) 

a.) is being used. The new 

(13) 



where the shorthand notation . . . ,a n ) 
POVM elements 

U n (a) = Gl(a)G n (a 
are normalized for all n: 



U(a.)dai ...a n = I , 



(14) 



as it follows from Eqs.(^Jl2.13). The probability of the whole sequence of 
outcomes can be written into the compact form: 



Pn{0-. 



tr 



n«(o-.)/5 



(15) 



The Eqs. (llj-lq) constitute a single (complicated) measurement. We 
invoke all considerations of state estimate from Sect. 2. In such a way shall 
we introduce the mixed state estimate 



Pn(°"-) 



n n (a.) 



(16) 



trII n (o-.) 

whose eigenstates, like in case of (^) , will be the pure state estimates. 
Same considerations that led to fidelities (||jl(]) in Sect. 2 apply invariably. 
We can, for instance, write the average fidelity in terms of the aposteriori 

J/2: 



state ( |ll|) emerging from a hypothetical apriori qubit state p- = p$ 

11 „ 

F n 



^ + ^Etr[^(a.)] 2 



(17) 



It is obvious that Fq = 1/2, and we expect F n is a monotone function of n. 
In Sect. 5 we prove that F n achieves the upper limit ([lj) even when each 
individual measurement is very unsharp. Actually, we shall prove that p' n 
tends to be pure for large n. The Sect. 4 prepares the mathematical tool of 
the proof. 
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4 Conditional master equation 



There is a particular class of sequential measurements which is treatable 
with good accuracy in terms of markovian stochastic differential equations. 
We assume long sequences of very unsharp measurements: 

n > 1 , A > 1 . (18) 



The asymptotic limit (T^, [10] 



n 

n, A — > co , = const (19) 

will be called the 'continuum limit'. In case of two-state systems, we assume 
that the measured observables &i,...,a n are Pauli-polarizations chosen in- 
dependently along random directions. Formally, let us count the succession 
of measurements as if they happened at constant rate v = 12/ A 2 . Accord- 
ingly, we replace the discrete parameter n by the continuous time: 

Yin . , 

*= ^2 ■ (2°) 

We consider all quantities as continuous functions of t, coarse-grained on 
scales S> \/v involving many measurements. In this limit an approximate 
theory emerges in the form of markovian stochastic differential equations. 
(The theory becomes exact in the continuum limit.) The aposteriori state, 



see Eq. (31) of the Appendix, satisfies the conditional (or selective) master 
equation: 

f = -^[^]] + {^-(^M , (21) 

where (a)t = tt[apt\. We have suppressed denoting the functional depen- 
dence of pt on the outcomes {<r T ;0 < r < t}. The wt is the standard 
isotropic white-noise and the equation must be interpreted in the sense of 
the Ito stochastic calculus. There is a second stochastic differential equation 
for the outcome: 

°t = 0)t + 2™t ■ (22) 
The features of the above equations have been well understood. In partic- 



ular, the solution p\ becomes asymptotically pure for long times [14, [lq| . 
This assures the saturation of average fidelity (|l~7| ) , as proven in the next 
section. So far, the stochastic differential equation governing the estimate 
p' t has been missing. We construct it in the Appendix. 
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5 Saturation of fidelity 



We are going to discuss the time dependence of the average fidelity Ft- Re- 
member that it corresponds to the (coarse-grained) n— dependent fidelity F n 
( |l7|) via t = 12ra/A 2 . The latter requires the knowledge of the hypothetical 
aposteriori state which, for a qubit, we shall parametrize by the polarization 
vector St = (<r) t : 

~? I + s t a 

Pt = 2 — ' ' Si ' = St ~ ' ( ' 

Recall that the initial state must be the hypothetical state 1/2 implying 
the initial value sq = 0. The stochastic 'master' equation (21) yields the 
following stochastic differential equation for the polarization vector: 

^ = -1st ~ 2 {s t w t ) s t + 2wt . (24) 
at 

This is an isotropic inhomogeneous spatial diffusion process. A stochastic 
differential equation for the squared norm (purity) follows from it: 

_L = 4(3 - s 2 t )(l - s 2 t ) + 4(1 - s 2 t )s t w t , (25) 

where wt is the standard white-noise. This is a one-dimensional inhomoge- 
neous diffusion. For long times the norm will approach the unity, therefore 
the aposteriori state becomes asymptotically pure. The author's Monte- 
Carlo calculations have shown that the purity s 2 is dominated by the drift 
term. Ignoring diffusion, the error remains within 2% and the analytic so- 
lution is possible: 

e 8t _ i 

E 4 = s 2 t = , t 1 /Q . (26) 



e 8t - 1/3 

Let us restore the original variable n = tA 2 /12 and substitute the above 
result into the expression ([D]) : 

11 11 e 96n/A2 - 1 

Fn= 2 + 6 E S * = 2 + 6e 9 WA2_ 1/3 " ^ 

The average fidelity approaches the optimum value 2 /3 after a characteristic 
number n ~ A 2 /96 of unsharp measurements. Recalling the conditions (|l8| ) 
we conclude that our result is valid for very unsharp measurements, i.e., A 
must be much greater than \/96 ~ 10. 
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6 Discussion 



We have discussed single quantum state determination via sequential non- 
projective (POVM) measurements in the limit of very unsharp measure- 
ments. We have proven that the known optimum average fidelity of esti- 
mating random qubits can be approached gradually with many successive 
random unsharp measurements. Whether this is true for non-random qubits 
is an open issue, but it is certainly tractable with the method of the present 
work. It may for instance turn out that one has to replace the strategy of 
random unsharp measurements by some adaptive strategy. 

We profited from analytic tools. We used the standard theory of (marko- 
vian) continuous quantum measurement and we completed it with the novel 
concept of continuous state estimation. The recent work p5| coined already 
a similar concept. It has, however, been restricted to the particular case of 
Gaussian states. Although we have detailed the concept for a single qubit, 
most of the equations are valid for any higher dimension N. The standard 
theory of continuous quantum measurement treats discrete and continuous 
observables on equal footing with the same formalism. We guess that also 
our continuous estimation formalism can be applied to the tomography of 



light quanta |16|], particularly to its Gaussian POVM formulation [17|. 

Stochastic differential equations, used so far for continuous measure- 
ment, will apply to optimum state determination as well. Continuous state 
determination is of interest every time when one is accumulating and ana- 
lyzing information from low rate quantum inference. These conditions are 
typical for an eavesdropper of secret quantum communication, a doner of 
n> 1 identical qubits into n + 1 identical qubits, or in tomography with 
low detection efficiency. 

I thank Nicolas Gisin for stimulating correspondence. This work was 
supported by the Hungarian OTKA Grant 32640. 



A Continuous measurement and estimation 



In the continuum limit (|l9|) , the outcome a. of sequential measurement 
(Sect. 3) becomes a (vectorial) function at of time. The basic mathemati- 
cal objects will be functionals of the outcome. First of all, we define the 
continuum limit of the sequential Kraus-operators (12) in terms of the 
time-ordered exponentials: 



G t [a] =Texp 



a T \ 2 dT 



(28) 



S 



The normalizing pre-factor of the exponential has been omitted and, as 
usual, will be incorporated in the functional measure d[a]. The above oper- 
ators yield the continuum limit of the sequential POVM (Il3|) : 



G\[a]G t [a] 



(29) 



It describes the isotropic continuous polarization measurement in the pe- 
riod [0, t]. The POVM satisfies the completeness relation at any time, with 
respect to the functional integration: 



/ 



u t [d)d[d] = i 



(30) 



The operators li|[cr], a kind of time-ordered Gaussian projectors, form a 
functional POVM for all time t. Given the random outcome {o>; < r < t}, 
the aposteriori state at time t takes this form: 



G t [d]pG\[a] 



ti- 



ll 



with the normalized functional probability distribution 



tr 



(31) 



(32) 



The Eqs. (|28|-[32]) constitute the model of isotropic continuous measurement 
of the polarization a. Similarly to the case of a single measurement, the 
choice of the estimate p' t is not unique. Following (0) and ([[(]) , as well 
as for mathematical convenience, we take 



Ptm 



trn t [a] 



(33) 



and, like in Sect. 2, we interprete it as the random average of its pure eigen- 
states. 

Contrary to the aposteriori state pt, the estimate state p' t does not sat- 
isfy an autonomous stochastic differential equation. Neither the composite 
object pt <8> p' t does. To construct a closed set of stochastic differential equa- 
tions, we introduce the state p' t where a hypothetical initial state 1/2 would 
have evolved to under the true operations Gt[cr] which the true apriori state 
Po = P had undergone: 

G t [a}G\[a} 



Pt 



tr fl t [a] 



(34) 
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Note in contrast to the preceeding sections, in particular to Sect. 2, that 
here we retain for p[ the probability (^) of the true continuous measure- 
ment. (Actually, we could have modified the notation p\.) We introduce 
two normalized variants of the Kraus-operators (28) : 



9t 



trU t p 



1/2 



9t 



1/2 



They will build up the time-dependent aposteriori (|3l|) 
, and the hypothetic state (p3) , respectively: 



(35) 

the estimate (|33|) 



Pt = gtpgi 



Pt 



tr/5 



= trpj 



ft 



(36) 



1 of these states follow from 



The normalizations trpt 
the normalizations (^) . Two time-dependent expectation values will be 
defined in function of the normalized operators ( |35| ) : 

(37) 
(38) 

For the sake of symmetry, I propose the normalized operators gt and g' t , 
yielding p t and p' t via (36) , to formulate a convenient couple of equations. 
An autonomous stochastic differential equation will exist for gt: 



0)t 


= tr 




( = tr 


Pt° 


0)1 


= ±tr Ugtftg't 


( = tr 


Pt° 



dg t 
dt 



1 



0)t\ 2 + 0-0)t)W t 



gt 



(39) 



This equation is equivalent with the well-known conditional master equation 
( f2l| ) . A new equation can be written down for g' t : 



dg\ 
dt 



\<7 



0)t 



(<?>I| 2 + I(<?h 



?|2 



+ (<? - 0)t)Wt 



g[ 



(40) 

This equation couples to the previous equation via [a)t in addition to the 
white- noise wt- The initial conditions are go = g'o = I- It is straightforward 
to show that the above equations preserve the normalization of pt and p' t . 

We outline the proof of the Eqs.( 39 ,[!(]). The proof will reside on the 
equation a t = (a)t + \wt of continuous measurement theory (22) . Let us 
substitute it into the definition (|28|) of the Kraus-operator G^[<t]. It yields 



Texp 



a — (a) T \ 2 dT +/(<?— (a) T )w T dr 



(41) 
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times a numeric factor which will be irrelevant for the normalized operators 
gt and g' t . It turns out that the above exponential is already the properly 
normalized gt- 



g t = T exp 



\a — {a) T \ 2 dr + 



(a) T )w T d,T 



(42) 



Indeed, differentiating the above equation yields exactly the Eq. 

Derivation of the novel equation (HO) for g' t is a bit more complicated. 
In addition to the exponential in Eq. (|42|) , we assume a further c-number 
differential for the sake of normalization (|35|) : 



a/ 

fJt 



exp 
Texp 



a T dr + T w r 



o 



t rt 
1, 



a — {a) T \ dr + (a — (a) T )w r dr 
n Jo 



(43) 



We calculate dg' t /dt and insert it into the normalization condition tr dp' t /dt = 0. 
This will yield the unique solutions 0t = {S)t — an( i a t = \Pt\ 2 - Inserting 
these results back into the equation of dg[/dt we obtain the Eq. ( |40| ) . 

The evolution of the aposteriori pt and the estimate state p' t is indirectly 
described by the coupled stochastic differential Eqs.(| 



4C). We mentioned 



that p obeys to a closed equation but pt®f>'t does not. From the above results 
it would be trivial to show that pt® g't® [g't) contains all information on 
pt ® p't and it does satisfy a closed stochastic differential equation. 
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